function f=dcdq2(c_jt,c,q_jt)
f=c_jt+c(1)*q_jt+c(2)*(q_jt.^2);
end
